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The QC Relaxation: A Theoretical and 
Computational Study on Optimal Power Flow 

Carleton Coffrin, Hassan L. Hijazi, and Pascal Van Hentenryck 


Abstract —Convex relaxations of the power flow equations and, 
in particular, the Semi-Definite Programming (SDP) and Second- 
Order Cone (SOC) relaxations, have attracted significant interest 
in recent years. The Quadratic Convex (QC) relaxation is a depar¬ 
ture from these relaxations in the sense that it imposes constraints 
to preserve stronger links between the voltage variables through 
convex envelopes of the polar representation. This paper is a 
systematic study of the QC relaxation for AC Optimal Power 
Flow with realistic side constraints. The main theoretical result 
shows that the QC relaxation is stronger than the SOC relaxation 
and neither dominates nor is dominated by the SDP relaxation. 
In addition, comprehensive computational results show that 
the QC relaxation may produce significant improvements in 
accuracy over the SOC relaxation at a reasonable computational 
cost, especially for networks with tight bounds on phase angle 
differences. The QC and SOC relaxations are also shown to be 
significantly faster and reliable compared to the SDP relaxation 
given the current state of the respective solvers. 

Index Terms —Optimization Methods, Convex Quadratic Op¬ 
timization, Optimal Power Flow 
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S 9 = p 9 + iq 9 
Co > Cl > C2 
»(•) 
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H 

x l ,x u 
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The set of nodes in the network 

The set of from edges in the network 

The set of to edges in the network 

imaginary number constant 

AC current 

AC power 

AC voltage 

Line impedance 

Line admittance 

Transformer properties 

Bus shunt admittance 

Product of two AC voltages 

Current magnitude squared, |/| 2 

Line charging 

Line apparent power thermal limit 
Phase angle difference limit 
AC power demand 
AC power generation 
Generation cost coefficients 
Real part of a complex number 
Imaginary part of a complex number 
Conjugate of a complex number 
Magnitude of a complex number, l 2 -norm 
Lower and upper bounds of x, respectively 
Convex envelope of x 
A constant value 
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I. Introduction 

C ONVEX relaxations of the power flow equations have 
attracted significant interest in recent years. They in¬ 
clude the Semi-Definite Programming (SDP) m, Second- 
Order Cone (SOC) (2), Convex-DistFlow (CDF) (3), and the 
recent Quadratic Convex (QC) E and Moment-Based a, 
a relaxations. Much of the excitement underlying this line 
of research comes from the fact that the SDP relaxation has 
shown to be tight on a variety of case studies |j7), opening 
a new avenue for accurate, reliable, and efficient solutions 
to a variety of power system applications. Indeed, industrial- 
strength optimization tools (e.g., Gurobi, cplex, Mosek) are 
now available to solve various classes of convex optimization 
problems. 

The relationships between the SDP, SOC, and CDF re¬ 
laxations is now largely well-understood: See 0, 0 for a 
comprehensive overview. In particular, the SOC and CDF 
relaxations are known to be equivalent and the SDP relaxation 
is at least as strong than both of these. However, little is 
known about the QC relaxation which is a significant departure 
from these more traditional relaxations. Indeed, one of the key 
features of the QC relaxation is to compute convex envelopes 
of the polar representation of the power flow equations in 
the hope of preserving stronger links between the voltage 
variables. This contrasts with the SDP and SOC relaxations 
which are derived from a lift-and-project approach on the 
complex representation. 

This paper fills this gap and provides a theoretical study of 
the QC relaxation as well as a comprehensive computational 
evaluation to compare the strengths and weaknesses of these 
relaxations. Our main contributions can be summarized as 
follows: 

1) The QC relaxation is stronger than the SOC relaxation. 

2) The QC relaxation neither dominates nor is dominated 
by the SDP relaxation. 

3) Computational results on optimal power flow show that 
the QC relaxation may bring significant benefits in 
accuracy over the SOC relaxation, especially for tight 
bounds on phase angle differences, for a reasonable loss 
in efficiency. 

4) The computational results also show that, with existing 
solvers, the SOC and QC relaxations are significantly 
faster and more reliable than the SDP relaxation. 

The theoretical results are derived using the equivalence of two 
classes of second-order cone constraints (in conjunction with 
the power equations), which provides an alternative formula¬ 
tion for the QC model which is interesting in its own right. 
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Moreover, to the best of our knowledge, the computational 
results also represent the most comprehensive comparison of 
these convex relaxations. They are obtained for optimal power 
flow problems with realistic side-constraints, featuring bus 
shunts, line charging, and transformers. 

The rest of the paper is organized as follows. Section [II] 
reviews the formulation of the AC-OPF problem from first 
principles and presents two equivalent formulations of this 
non-convex optimization problem. Section [HI] derives the SDP, 
QC, and SOC relaxations. Section |rv| illustrates their behavior 
on a well-known 3-bus example. Section [V] presents an alter¬ 
native formulation of the QC relaxation which is a convenient 
tool for subsequent proofs. Section [VI] presents the theoretical 
results linking the QC to the other relaxations. Section |VII| 
reports the computational results for the three relaxations on 
93 AC-OPF test cases, and Section VIII concludes the paper. 


II. AC Optimal Power Flow 

This section reviews the specification of AC Optimal Power 
Flow (AC-OPF) and introduces the notations used in the paper. 
In the equations, constants are always in bold face. The AC 
power flow equations are based on complex quantities for 
current I, voltage V , admittance Y, and power S, which are 
linked by the physical properties of Kirchhoff’s Current Law 
(KCL), i.e., 


I?-It= E I l3 (1) 

( i,j)eEuE R 

Ohm’s Law, i.e., 

u 3 = Y l] (y l -v 3 ) (2) 

and the definition of AC power, i.e., 

Si, = Vil* 3 (3) 

Combining these three properties yields the AC Power Flow 
equations, i.e., 

S? -Sf= E S U yieN (4a > 

( i,j)eEuE R 

S l3 =Y*V Z V* -Y*V t V* (■ i,j)eEuE R (4b) 

These non-convex nonlinear equations define how power flows 
in the network and are a core building block in many power 
system applications. However, practical applications typically 
include various operational side constraints on the power flow. 
We now review some of the most significant ones. 

Generator Capabilities: AC generators have limitations 
on the amount of active and reactive power they can produce 
S 9 , which is characterized by a generation capability curve 
ED- Such curves typically define nonlinear convex regions 
which are typically approximated by boxes in AC transmission 
system test cases, i.e., 

sf < s? Sf u Vi e N (5a) 


Line Thermal Limit: AC power lines have thermal limits 
m to prevent lines from sagging and automatic protection 
devices from activating. These limits are typically given in 
Volt Amp units and constrain the apparent power flows on the 
lines, i.e., 

V(i,j)eEuE R ( 6 ) 

Bus Voltage Limits: Voltages in AC power systems should 
not vary too far (typically ±10%) from some nominal base 
value ED- This is accomplished by putting bounds on the 
voltage magnitudes, i.e., 

vl ^ \Vi\ ^ Vi e N (7) 

A variety of power flow formulations only have variables 
for the square of the voltage magnitude, i.e., |V)| 2 . In such 
cases, the voltage bound constrains can be incorporated via 
the following constraints: 

(v‘) 2 ^ M 2 < («?) 2 VieN (8) 

Phase Angle Differences: Small phase angle differences 
are also a design imperative in AC power systems Ql and it 
has been suggested that phase angle differences are typically 
less than 10 degrees in practice HD- These constraints have 
not typically been incorporated in AC transmission test cases 
b a. However, recent work na, a have observed that 
incorporating Phase Angle Difference (PAD) constraints, i.e., 

(9) 

is useful in the convexification of the AC power flow equa¬ 
tions. For simplicity, this paper assumes that the phase an¬ 
gle difference bounds are symmetrical and within the range 
(—7 t/2, 7t/ 2), i.e., 

(i,j)eE (10) 

but the results presented here can be extended to more general 
cases. Observe also that the PAD constraints © can be 
implemented as a linear relation of the real and imaginary 
components of ViV* HU, i.e. V(i,j) e E, 

tan(~0$)X (ViV*) (VS^*) *St«m(0£)K (W/) (11) 

The usefulness of this formulation will be apparent later in the 
paper. 

Other Constraints: Other line flow constraints have been 
proposed, such as, active power limits and voltage difference 
limits 13. El. However, we do not consider them here since, 
to the best of our knowledge, test cases incorporating these 
constraints are not readily available. 

Objective Functions: The last component in formulating 
OPF problems is an objective function. The two classic 
objective functions are line loss minimization, i.e., 

minimize: E 9ft(S?) (12) 

ieN 

and generator fuel cost minimization, i.e., 

minimize: E c 2 i(9ft(S?)) 2 + d MS?) + c 0l (13) 

ieN 

Observe that objective © is a special case of objective © 
where C 2 i = 0, cu = 1, Co,: = 0 ( ieN ) fl5l . Hence, the rest of 
this paper focuses on objective ([T3|. 
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Model 1 AC-OPF 

variables: Sf (Vi e N), VfiVi e N) 

minimize: c 2i ()R.(Sf)) 2 + c-|,; s ft(,S'f) + c 0i (15a) 

ieN 

subject to: 

< |VJ] < < Vi e N (15b) 

sf < Sf Sf u Vi eN (15c) 

|Sy|<«J5 V(i,j)eEuE R (15d) 

Sf - Sf = ^ Sij VieN (15e) 

( i,j)eEuE R 

Sij = Y*-ViV* — Yi*ViV* (i,j)eEvE R (15f) 
- 0% < AVIV?) ^ 0$ V(i,j) g E (15g) 


Model 3 The SDP Relaxation AC-OPF-W-SDP. 

variables: S'?(Vi g N), Wij(Vi,j g N) 

minimize: ( |16a| ) 
subject to: ( | 16c| )— ( fl~6l] ) 

W > 0 (17a) 


Model 4 The SOC Relaxation AC-OPF-W-SOC. 

variables: Sf(V« g N), Wij(V(i,j) e E), Wu(Vi g N) : real 

minimize: ( |16a| ) 
subject to: ( | 16c| )— ( p~6T| ) 

\W tJ \ 2 ^ WaWjj V(i,j)eE (18a) 


Model 2 AC-OPF-W 


variables: S?(VieN), Vi(VieN), Wij(Vi,j 

g N) 

minimize: ^ c 2 i(3?(S'f)) 2 + ci + c 0i 

(16a) 

isN 


subject to: 


Wij = ViV? Vi g N, Vj g N 

(16b) 

(it-) 2 < Wu < «) 2 Vi e N 

(16c) 

Sf 1 S? Sf u VieN 

(16d) 

Sf - sf = ^ Sij Vi g N 

(16e) 

(i,j)eEuE R 


S ij = Y*W ii -Y* j W ij (i,j) g E 

(16f) 

Sji = Y*jWjj — Y*jW? (Lj) G /•; 

(16g) 

\Sij\ < K) V(i,j) e E v E r 

(16h) 


tan (-G^AWij) ^(Wy) < tan(0£)* l(W zj ) (16i) 

V(*, j) G E 


AC-OPF: Combining the AC power flow equations, the 
side constraints, and the objective function, yields the well- 
known AC-OPF formulation presented in Model |T| Observe 
that, in Model [T] the non-convexities arises solely from the 
product of the voltages (i.e., V, V*) and they can be isolated 
by introducing new W variables to represent the products of 

vs mi, 0, (ra, m, i.e, 

ViV* = Wij (i,j g N). (14) 

Model [2] presents an equivalent version of the AC-OPF, where 
the W factorization has been incorporated and the only source 
of non-convexity is in constraint ( | 16b| ). Note that this section 
has introduced the simplest form of the AC-OPF problem and 
that real-world applications feature a variety of extensions 
as discussed at length in Q3, Go). In practice, this non- 
convex nonlinear optimization problem is typically solved with 
numerical methods (e.g. IPM, SLP) GTl , 1221 . which provide 
locally optimal solutions if they converge to a feasible point. 


III. Convex Relaxations of Optimal Power Flow 


Since the AC-OPF problem is NP-Hard lf23l . Il24l and 
numerical methods provide limited guarantees for determining 
feasibility and global optimally, significant attention has been 
devoted to finding convex relaxations of Model |T| Such 
relaxations are appealing because they are computationally 
efficient and may be used to: 


1) bound the quality of AC-OPF solutions produced by 
locally optimal methods; 

2) prove that a particular AC-OPF problem has no solution; 

3) produce a solution that is feasible in the original non- 
convex problem Q, thus solving the AC-OPF and 
guaranteeing that the solution is globally optimal. 


The ability to provide bounds is particularly important for the 
numerous mixed-integer nonlinear optimization problems that 
arise in power system applications. For these reasons, a variety 
of convex relaxations of the AC-OPF have been developed in¬ 
cluding, the SDP [[[), QC 0|, SOC GL and Convex-DistFlow 
0, which are reviewed in detail in this section. Moreover, 
since the SOC and Convex-DistFlow relaxations have been 
shown to be equivalent 11251 . this paper focuses on the SDP, 
SOC, and QC relaxations only and shows how they are derived 
from Model [2] The key insight is that each relaxation presents 

1 6bb, which 


a different approach to convexifing constraints 
are the only source of non-convexity in Model [2 

The Semi-Definite Programming (SDP) Relaxation: ex¬ 
ploits the fact that the W variables are defined by V(V*) T , 
which ensures that W is positive semi-definite (denoted by 
W > 0) and has rank i CD, 0, OH. These conditions are 
sufficient to enforce constraints ( fl6b| 1261 . i.e., 


Wij = VW* (i, j g N) <=> W > 0 A rank(FF) = l 

The SDP relaxation G3. GU then drops the rank constraint 
to obtain Model [3] 

The Second Order Cone (SOC) Relaxation: convexifies 
each constraint of ( |16b| ) separately, instead of considering them 
globally as in the SDP relaxation. The SOC relaxation takes 
the absolute square of each constraint, refactors it, and then 
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relaxes the equality into an inequality, i.e.. 


Wu = v,v; 

(19a) 

WijW* = ViV*V*V 3 

(19b) 

IWyf = WjjW^ 

(19c) 

\Wi il 2 < WjjWtf 

(19d) 


Equation ( |19d| > is a rotated second-order cone constraint which 
is widely supported by industrial optimization tools. It can, in 
fact, be rewritten in the standard form of a second-order cone 
constraint as. 


( 2 W t , \ 

\W« ~ W n) 


< Wu + Wjj 


( 20 ) 


The complete SOC formulation is presented in Model [4] Note 
that this relaxation requires fewer W variables than Model [3] 
Due to the sparsity of AC power networks, this size reduction 
can lead to significant memory and computational savings. 

The Quadratic Convex (QC) Relaxation: was introduced 
to preserve stronger links between the voltage variables Bl¬ 
it represents the voltages in polar form (i.e., V = vZ.6) and 
links these real variables to the W variables, along the lines 
of Q6), (nl, ESI, ED, using the following equations: 


Wu = vf i e N (21a) 

ft(Wy) = ViVj cos(6»j - 6j) V(i,j) g E (21b) 

sin(0j - 9j) g E (21c) 


The QC relaxation then relaxes these equations by taking 
tight convex envelopes of their nonlinear terms, exploiting the 
operational limits for v t , v 3 , 0, — 0 t . The convex envelopes for 
the square and product of variables are well-known (30l . i.e.. 


(xV = 


X ^ X 2 

X < ( x u + x l )x — x u x l 


(xy) M - 


'xy ^ x l y + y l x — x l y l 
xy > x u y + y u x — x u y u 
xy ^ x y + y u x — x y u 
Jicy x u y + y l x — x u y l 


(T-CONV) 


(M-CONV) 


Under our assumptions that the phase angle bound satisfies 
0 ^ 0 A s; ^ and is symmetric, convex envelopes for sine 
(S-CONV) and cosine (C-CONV) (4| are given by. 


Ux < cos (sf ) (x - + sin 

[SS ^ cos (2f ) (x + - sin 


<COs(a:)) L = 


r -1 1—cos(a3 lt ) 2 
CX < 1 - X 

cx ^ cos(x u ) 


In the following, we abuse notation and also use (/(-)) C 
to denote the variable on the left-hand side of the convex 
envelope C for function /(•). When such an expression is used 
inside an equation, the constraints (/(-)) C are also added to 
the model. 

Convex envelopes for equations \2 1 a| >—( |2Tc] > can be obtained 
by composing the convex envelopes of the functions for 


Model 5 The QC Relaxation AC-OPF-C-QC. 

variables: 5f(V* e N), Wij(V(i,j) e E), Wjj (Vi e N) : real 
Vi Z0 4 (Vi g N), kMi,j) g E) 

minimize: (16a| ) 
subject to: (16c| )-( fT6il > 


II 

ro 

! > T ieN 



(22a) 

w = 

({ViVj) M <cos(6'j - 

e 3 )) c ) M 

v(f,i) g E 

(22b) 

3(W y ) = 

((viVj) M (sm(6i - 


V(*,i)G£ 

(22c) 

Sij + Sji 



(22d) 

ISjj| 2 < WjjZjj V(i, j) g E 



(22e) 



Fig. 1. A Network Diagram for the 3-Bus Example. 


square, sine, cosine, and the product of two variables, i.e., 

Wj = (v 2 ) T i g N (23a) 

cos(0j - 0,)> C > M V(i, j) G E (23b) 
= «^^> M <sin(^j - %)> S > M V(i, j) g E (23c) 

The QC relaxation also proposes to strengthen these convex 
envelopes with a second-order cone constraint based on the 
absolute square of line power flow (3]l, first proposed in 0. 
This requires a new variable U 3 for each line (i,j) G E 
that captures the current magnitude squared on that line. The 
following constraints are added to link the Uj variables to the 
existing model variables. 

S%j + Sji = Zijly V(i,j) G E (24a) 

|Sy| 2 < Wijiy V(i,j)eE (24b) 

The complete QC relaxation is presented in Model [5] This 
model is annotated as C-QC, as the second-order cone con¬ 
straints use current variables. The motivation for this distinc¬ 
tion will become clear in Section [V] 

IV. An Illustrative Example 
This section illustrates the three main power flow relaxations 
on the 3-bus network from ED, which has proven to be an 
excellent test case for power flow relaxations. This system is 
depicted in Figure [T] and the associated network parameters are 
given in Table [I] This network is designed to have very few 
binding constraints. Hence, the generator and line limits are 
set to large non-binding values, except for the thermal limit 
constraint on the line between buses 2 and 3, which is set to 50 
MVA. In addition to its base configuration, we also consider 
this network with reduced phase angle difference bounds of 
18°. IPOPT (32l is used as a heuristic 1331 to find a feasible 
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TABLE I 

Three-Bus System Network Data (100 MVA Base). 


Bus Parameters 

Bus 

P d 

qd 

v l 

v u 

1 

no 

40 

0.9 

1.1 

2 

no 

40 

0.9 

1.1 

3 

95 

50 

0.9 

1.1 


Line Parameters 

From-To Bus 

r 

X 

b c 

s u 

0* 

1-2 

0.042 

0.90 

0.30 

00 

30° 

2-3 

0.025 

0.75 

0.70 

50 

o 

o 

CO 

1-3 

0.065 

0.62 

0.45 

00 

30° 


Generator Parameters 

Generator 

p 9« 

qQi ^ q9 u 

C 2 

Cl 

CO 

i 

0, oo 

—oo, oo 

0.110 

5.0 

0 

2 

0, oo 

—00, 00 

0.085 

1.2 

0 

3 

0,0 

—oo, oo 

0 

0 

0 


TABLE II 

AC-OPF Bounds using Relaxations on the 3-Bus Case. 


Test Case 

$/h 

AC 

Optim 

SDP 

ality Gc 
QC 

ip (%) 
SOC 

Base 

5812 

0.39 

1.24 

1.32 

6» A = 18° 

5992 

2.06 

1.24 

4.28 


solution to the AC-OPF and we measure the optimally gap 
between the heuristic and a relaxation using the formula 
Heuristic — Relaxation 
Heuristic 

Table |h| summarizes the results^] In the base configuration, the 
SDP relaxation has the smallest optimality gap. In the d A = 
18° case, the QC relaxation has the smallest optimality gap, 
while reducing the bound on phase angle differences increases 
the optimality gap for both the SDP and SOC relaxations. 
This small network highlights two important results. First, the 
SDP relaxation does not dominate the QC relaxation and vice- 
versa. Second, the SDP and QC relaxations dominate the SOC 
relaxation. The next two sections prove that this last result 
holds for all networks. 

V. An Alternate Form of the QC Relaxation 

Section m introduced two types of second-order cone 
constraints. Model [4] uses a SOC constraint based on the 
absolute square of the voltage product 12, i.e., 

| Wijf^WiiWjj (25) 

while Model 0 uses a SOC constraint based on the absolute 

square of the power flow 0, i.e., 

\Sijl 2 < Wukj. (26) 

We now show that, in conjunction with the power flow equa¬ 
tions ([T6fHT6B, these two SOC formulations are equivalent. 
More precisely, we show that 

Sa = Y*W« - Y*Wij (i,j) £ E 

Sp = Y*W 30 -Y*W* (i,j) e E (W-SOC) 

m 3 \ 2 WaWjj (i,j) £ E 


1 On this small example a nonlinear global optimization solver was used to 
prove that the heuristic solutions are in fact globally optimal. Such a validation 
is not possible on larger test cases. 


Model 6 The Alternate QC Relaxation AC-OPF-W-QC 

variables: S?(V* e N), Wij(V(i,j) e E), Wu(\/i e N) : real 
UiZ0i(Vi g N ) 

minimize: d 16a| ) 

subject to: ( |16c| )-( [T6i] l, ( |22a[ i-( f22c| i. ( |18a| ) 


is equivalent to 

Sij = Y*w ii -Y*w ij (i.j) C /•: 
Sa = Y£W jj -Y*WZ (i.j) e /•; 

&ij T ^ij Ljj ( i , j ) G E 

I Sij\ 2 (i,j) g E. 


(C-SOC) 


This equivalence suggests an alternative formulation of the 
QC relaxation which is given in Model [6] and establishes a 
clear connection between Models [4] and [5] Throughout this 
paper, we use W and C to denote which of these equivalent 
formulations is used. 

We now prove these results. The following lemma, whose 
proof is straight-forward and can be found in the Appendix, 
establishes some useful equalities. 


Lemma V.l. The following four equalities hold: 

1 ) \S tJ \ 2 = \Yij\ 2 (Wl - WaWij - W vl W* + \Wij\ 2 ). 

2) \Wij\ 2 = Wl - W*Z*S i3 Wi.ZijSfj + \Z tj \ 2 1%| 2 . 

3) kj = \Y VJ \ 2 (W U + W 3j - U'„ - W*). 

4) W j3 = W u - Z*S l3 - Z l3 S* + |%|V 

We are now ready to prove the main result of this section. 
Theorem V.2. ( |C-SQC[ i is equivalent to ( |W-SOC| ). 


Proof: 

The proof is similar in spirit to those presented in 

C-SOC ' Every solution to ( W-SOC| l is a 


W-SOC 


solution to ( C-SOC| l. Given a solution to ( |W-SOQ , by equality 
(3) in Lemma V.l we assign lij as follows: 


kj = \Yj\ 2 (w u - w i3 - w* + w jd ) (i,j) g /•: 


This assignment satisfies the power loss constraint (J22dj) by 
definition of the power. It remains to show that second-order 
cone constraint in ( C-SOQ is satisfied. Using equalities (1) 
and (3) in Lemma [Vi] we obtain 


l^f = l^il 2 {Wl - WaWij - W vl W* + \W t f 2 ) 

I Sij \ 2 Z \Y i3 \ 2 (Wl - WaWij - W U W* + WaWjj) 
\Sij\ 2 < WalYf 2 (IF;,- - Wij - IF*- + IF;;) 

\Sij\ 2 < Wakj. 


|C-5DC| W-SOC\ Every solution to ( |C-SOC[ i 


is a 


solution to ( |W-SOC >. We show that the values of W KJ in 
( |C-SOC| > satisfy the second-order cone constraint in ( |W-SOC| . 
Using equalities (2) and (4) in Lemma V.l and the fact that 
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Fig. 2. A Venn Diagram of the Solutions Sets for Various AC Power Flow 
Relaxations (set sizes in this illustration are not to scale). 


Wu = W? since Wa is a real number, we have 

\Wij \ 2 = wl - W*zf 3 s l3 - WuZiiSti + 

Wf 2 < w 2 - W*Z*S l3 - WiiZijS* + z„ Hv n l,j 
\W l3 \ 2 < WuiWu - Z*S,j - Z,jSf :i + \Zij\%j) 

\W tJ \ 2 < WiiWjj 

and the result follows. ■ 


Corollary V.3. Model [5] is equivalent to Model [6] 

Computational results on these two formulations are presented 
in the Appendix. The main message is that the |C-SOC| 
formulation is preferable to W-SOC in the current state of 
the solving technology, especially on very large networks. 

It is important to note that, for clarity, the proofs are pre¬ 
sented on the purest version of the AC power flow equations. 
Transmission system test cases typically include additional 
parameters such as bus shunts, line charging, and transformers. 
Proofs that these results can be extended to include the 
additional parameters in transmission system test cases are 
presented in the Appendix. 


VI. Relations of the Power Flow Relaxations 

We are now in a position to state the relationships between 
the convex relaxations. Recall that model M\ is a relaxation of 
model M 2 , denoted by M 2 £ Mi, if the solution set of M 2 is 
included in the solution set Mi. We use Mi A M 2 to denote 
the fact that neither M 2 £ Mi nor Mi c M 2 holds. Since 
our relaxations have different sets of variables, we define the 
solution set as the assignments to the Wij variables. 

Theorem VI.l. The following properties, illustrated in Figure 
[2] hold: 

1) SDP c SOC. 

2) SDP A QC. 

3) QC c SOC. 

Proof: Properties (1) and (2) follows from m and 
Section [Tv| respectively. For Property (3), observe that the set 
of constraints in Model [6] (W-QC) is a superset of those in 
Model [4] The result follows from Corollary |V.3| ■ 

Observe that the additional constraints ( |22a| i-( |22cj ) in the QC 
formulations are parameterized by the 0 A . As 0 A grows larger, 
the QC model reduces to the SOC model. Clearly, the strength 
of the QC relaxation is sensitive to this input parameter, as 
illustrated in Section IfVl 


VII. Computational Evaluation 

This section presents a computational evaluation of the 
relaxations and address the following questions: 

1) How big are the optimality gaps in practice? 

2) What are the runtime requirements of the relaxations? 

3) How robust is the solving technology for the relaxations? 
The relaxations were compared on 105 state-of-the-art AC- 
OPF transmission system test cases from the NESTA v0.4.0 
archive (35). These test cases range from as few as 3 buses to 
as many as 9000 and consist of 35 different networks under a 
typical operating condition (TYP), a congested operating con¬ 
dition (API), and a small angle difference condition (S AD)0 

Experimental Setting: All of the computations are con¬ 
ducted on Dell PowerEdge R415 servers with Dual 2.8GHz 
AMD 6-Core Opteron 4184 CPUs and 64GB of memory. 
IPOPT 3.12 (32) with linear solver ma27 (36l , as suggested by 
EH, was used as a heuristic for finding locally optimal feasi¬ 
ble solutions to the non-convex AC-OPF formulated in AMPL 
[|38ll . The SDP relaxation was executed on the state-of-the-art 
implementation (39) which uses a branch decomposition ]40l 
with a minor extension to add constraint ( |16i| >. The SDP solver 
SDPT3 4.0 SD was used with the modifications suggested 
in (39l . The second-order cone models were formulated in 
AMPL and IPOPT was used to solve the models. Numerical 
stability appears to be a significant challenge on the power 
networks with more than 1000 buses (42). Note that IPOPT 
is single-threaded and does not take advantage of the multiple 
cores available in the computation servers. This gives some 
computational advantage to the SDP solver, which utilizes 
multiple cores. 

Challenging Test Cases: We observe that 52 of the 105 
test cases considered have an optimality gap of less than 1.0% 
with the SOC relaxation. Such test cases are not particularly 
useful for this study as the improvements of the SDP and QC 
models are minor. Hence, we focus our attention on the 53 
test cases where the SOC optimality gap is greater than 1.0%. 
The results are displayed in Table [III] 


A. The Quality of the Relaxations 

The first six columns of Table III present the optimality 
gaps for each of the relaxations on the 53 challenging NESTA 
test cases. Note that bold values indicate cases where the 
relaxation produced a solution to the non-convex AC power 
flow problem. The table illustrates that this is a rare occurrence 
in the cases considered. 

The SDP Relaxation: Overall, the SDP relaxation tends 
to be the tightest, often featuring optimality gaps below 1.0%. 
In 5 of the 53 cases, the SDP relaxation even produces a 
feasible AC power flow solution (as first observed in (71). 
However, with six notable cases where the gap is above 5%, 
it is clear that small gaps are not guaranteed. In some cases, 
the optimality gap can be as large as 30%. 

A significant issue with the SDP relaxation is the reliability 
of the solving technology. Even after applying the solver 


-Nine test cases based on the EIR Grid network were omitted from 
evaluation because the AC-OPF-W-SDP solver did not support inactive buses. 
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TABLE III 

Quality and Runtime Results of AC Power Flow Relaxations 


Test Case 

$/h 

AC 

Optimality Gap (%) 

SDP | QC | SOC | CP 

Runtime (seconds) 

AC | SDP | QC | SOC | CP 

Typical Operating Conditions (TYP) 

nesta case3 lmbd 

5812.64 

0.39 

1.24 

1.32 

2.99 

0.12 

4.16 

0.07 

0.05 

0.03 

nesta case5 pjm 

17551.89 

5.22 

14.54 

14.54 

15.62 

0.04 

5.36 

0.09 

0.03 

0.05 

nesta case30 ieee 

204.97 

0.00 

15.64 

15.88 

27.91 

0.09 

8.38 

0.17 

0.07 

0.06 

nesta casel 18 ieee 

3718.64 

0.06 

1.72 

2.07 

7.87 

0.41 

12.62 

0.87 

0.43 

0.05 

nesta case 162 ieee dtc 

4230.23 

1.08 

4.00 

4.03 

15.44 

0.61 

35.20 

1.48 

0.31 

0.04 

nesta case300 ieee 

16891.28 

0.08 

1.17 

1.18 

n.a. 

0.80 

29.69 

2.83 

0.65 

n.a. 

nesta case2224 edin 

38127.69 

1.22 

6.03 

6.09 

8.45 

11.42 

690.16 

65.59 

45.99 

0.33 

nesta case2383wp mp 

1868511.78 

0.37 

1.04 

1.05 

5.35 

12.41 

1966.10 

57.87 

12.91 

0.80 

nesta case3012wp mp 

2600842.72 

— 

1.00 

1.02 

n.a. 

12.40 

14588.79t 

53.59 

19.15 

n.a. 

nesta case9241 pegase 

315913.26 

— 

1.67 

— 

n.a. 

132.25 

— 

3064.42 

— 

n.a. 

Congested Operating Conditions (API) 

nesta case3 lmbd api 

367.74 

1.26 

1.83 

3.30 

14.79 

0.18 

4.41 

0.09 

0.05 

0.23 

nesta case6 ww api 

273.76 

0.00* 

13.14 

13.33 

17.17 

0.34 

13.19 

0.07 

0.06 

0.03 

nesta case 14 ieee api 

325.56 

0.00 

1.34 

1.34 

8.89 

0.19 

5.64 

0.11 

0.08 

0.94 

nesta case24 ieee rts api 

6421.37 

1.45 

13.77 

20.70 

24.12 

0.14 

7.50 

0.26 

0.09 

0.04 

nesta case30 as api 

571.13 

0.00 

4.76 

4.76 

8.01 

0.38 

6.12 

0.17 

0.11 

1.11 

nesta case30 fsr api 

372.14 

11.06 

45.97 

45.97 

48.80 

0.25 

7.25 

0.19 

0.09 

0.92 

nesta case30 ieee api 

415.53 

0.00 

1.01 

1.01 

12.75 

0.07 

6.60 

0.19 

0.09 

0.03 

nesta case39 epri api 

7466.25 

0.00 

2.97 

2.99 

13.31 

0.10 

7.36 

0.29 

0.12 

0.04 

nesta case73 ieee rts api 

20123.98 

4.29 

12.01 

14.34 

17.83 

0.48 

10.03 

0.66 

0.20 

0.06 

nesta case89 pegase api 

4288.02 

18.11 

20.39 

20.43 

22.60 

1.16 

21.58 

1.29 

0.81 

0.04 

nesta case 118 ieee api 

10325.27 

31.50 

43.93 

44.08 

49.69 

0.46 

12.59 

0.84 

0.25 

0.05 

nesta case 162 ieee dtc api 

6111.68 

0.85 

1.33 

1.34 

19.39 

0.50 

36.85 

1.53 

0.39 

0.05 

nesta case 189 edin api 

1982.82 

0.05 

5.78 

5.78 

n.a. 

1.07 

16.10 

1.14 

0.33 

n.a. 

nesta case2224 edin api 

46235.43 

1.10 

2.77 

2.77 

9.07 

12.28 

672.04 

81.66 

88.33 

0.33 

nesta case2383wp mp api 

23499.48 

0.10 

1.12 

1.12 

3.10 

9.50 

1421.39 

28.37 

10.25 

0.34 

nesta case2736sp mp api 

25437.70 

0.07 

1.32 

1.33 

3.89 

9.21 

2278.77 

41.29 

10.51 

0.36 

nesta case2737sop mp api 

21192.40 

0.00 

1.05 

1.06 

4.62 

9.29 

1887.22 

30.94 

9.91 

0.32 

nesta case2869 pegase api 

96573.10 

0.92* 

1.49 

1.49 

5.16 

21.03 

1579.87 

102.55 

161.96 

0.37 

nesta case3120sp mp api 

22874.98 

— 

3.02 

3.03 

n.a. 

14.92 

15018.93T 

41.72 

12.19 

n.a. 

nesta case9241 pegase api 

241975.18 

— 

2.45 

2.59 

n.a. 

140.73 

— 

3511.60 

8387.11 

n.a. 

Small Angle Difference Conditions (SAD) 

nesta case3 lmbd sad 

5992.72 

2.06 

1.24* 

4.28 

5.90 

0.19 

4.39 

0.10 

0.05 

0.03 

nesta case4 gs sad 

324.02 

0.05 

0.81 

4.90 

66.06 

0.24 

4.16 

0.06 

0.06 

0.07 

nesta case5 pjm sad 

26423.32 

0.00 

1.10 

3.61 

43.95 

0.08 

5.35 

0.11 

0.05 

0.03 

nesta case6 c sad 

24.43 

0.00 

0.40 

1.36 

6.79 

0.26 

5.32 

0.11 

0.05 

0.02 

nesta case9 wscc sad 

5590.09 

0.00 

0.41 

1.50 

6.69 

0.14 

4.18 

0.19 

0.05 

0.03 

nesta case24 ieee rts sad 

79804.96 

6.05 

3.88 

11.42 

23.56 

0.10 

6.24 

0.30 

0.11 

0.04 

nesta case29 edin sad 

46933.26 

28.44 

20.57 

34.47 

36.79 

0.70 

9.19 

1.73 

0.27 

0.06 

nesta case30 as sad 

914.44 

0.47 

3.07 

9.16 

16.06 

0.18 

6.49 

0.22 

0.09 

0.03 

nesta case30 ieee sad 

205.11 

0.00 

3.96 

5.84 

27.96 

0.12 

7.49 

0.18 

0.09 

0.03 

nesta case7 3 ieee rts sad 

235241.70 

4.10 

3.51 

8.37 

22.21 

0.30 

9.48 

0.87 

0.20 

0.07 

nesta casel 18 ieee sad 

4324.17 

7.57 

8.32 

12.89 

20.77 

0.56 

14.14 

0.98 

0.31 

0.06 

nesta case 162 ieee dtc sad 

4369.19 

3.65 

6.91 

7.08 

18.13 

0.81 

39.71 

1.70 

0.36 

0.05 

nesta case 189 edin sad 

914.61 

1.20* 

2.22 

2.25 

n.a. 

0.65 

14.83 

1.27 

0.46 

n.a. 

nesta case300 ieee sad 

16910.23 

0.13 

1.16 

1.26 

n.a. 

1.01 

29.63 

2.81 

0.76 

n.a. 

nesta case2224 edin sad 

38385.14 

1.22 

5.57 

6.18 

9.06 

11.53 

691.53 

50.34 

65.68 

0.33 

nesta case2383wp mp sad 

1935308.12 

1.30 

2.97 

4.00 

8.62 

16.25 

1785.26 

40.71 

12.57 

0.80 

nesta case2736sp mp sad 

1337042.77 

2.18* 

2.01 

2.34 

4.56 

13.22 

1737.25 

35.42 

11.31 

0.48 

nesta case2737sop mp sad 

795429.36 

2.24* 

2.21 

2.42 

3.95 

13.01 

2153.37 

32.05 

9.69 

0.39 

nesta case2746wp mp sad 

1672150.46 

2.41* 

1.83 

2.44 

5.43 

14.01 

2840.32 

35.66 

13.32 

0.56 

nesta case2746wop mp sad 

1241955.30 

2.71* 

2.48 

2.94 

5.14 

14.51 

2306.18 

32.41 

23.22 

0.42 

nesta case3012wp mp sad 

2635451.29 

— 

1.92 

2.12 

n.a. 

15.79 

13548.13^ 

46.59 

28.41 

n.a. 

nesta case3120sp mp sad 

2203807.23 

— 

2.56 

2.79 

n.a. 

30.01 

16804.55T 

53.81 

15.69 

n.a. 

nesta_case9241_pegase_sad 

315932.06 

— 

0.80 

1.75 

n.a. 

80.30 

— 

3531.62 

33437.86 

n.a. 


bold - the relaxation provided a feasible AC power flow, * - solver reported numerical accuracy warnings, —,f - iteration or memory limit 


modifications suggested in 1391 , the solver fails to converge to 
a solution before hitting the default iteration limit on 8 of the 
53 test cases shown, it reports numerical accuracy warnings 
on 4 of the test cases, and ran out of memory on the 3 test 
cases with more 9000 nodes. 

The QC and SOC Relaxations: As suggested by 
the theoretical study in Section VI when the phase an¬ 
gle difference bounds are large, the QC relaxation is 


quite similar to the SOC relaxation. However, when the 
phase angle difference bounds are tight (e.g., in the 
SAD cases), the QC relaxation has significant benefits 
over the SOC relaxation. On average, the SDP relaxation 
dominates the QC and SOC relaxations. However, there 

are several notable cases (e.g. nesta_case24_ieee_rts_sad, 

nesta_case29_edin_sad, nesta_case73_ieee_rts_sad) where 
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the QC relaxation dominates the SDP relaxation. 

The Copper Plate (CP) Relaxation: This relaxation in¬ 
dicates the cost of supplying power to the loads when there 
are no line losses or network constraints (43) . and is included 
in the table as a point of reference. Note that this relaxation 
cannot be applied to networks containing lines with negative 
resistance or impedance, as indicated by “n.a.”. 


B. The Performance of the Relaxations 

Detailed runtime results for the heuristic solution method 
and the relaxations are presented in the last four columns of 
Table III The AC heuristic is fast, often taking less than 1 
second on test cases with less than 1000 buses. The SOC 
relaxation most often has very similar performance to the AC 
heuristic. The additional constraints in the QC relaxation add 
a factor 2-5 on top of the SOC relaxation. In contrast to these 
other methods, the SDP relaxation stands out, taking 10-100 
times longer. It is interesting to observe, in the 5 cases where 
the SDP relaxation finds an AC-feasible solution, the heuristic 
finds a solution of equal quality in a fraction of the time. 
Focusing on the test cases where the SDP fails to converge, 
we observe that the failure occurs after several minutes of 
computation, further emphasizing the reliability issue. 


VIII. Conclusion 

This paper compared the QC relaxation of the power flow 
equations with the well-understood SDP and SOC relaxations 
both theoretically and experimentally. Its two main contribu¬ 
tions are as follows: 

1) The QC relaxation is stronger than the SOC relaxation 
and neither dominates nor is dominated by the SDP 
relaxation. 

2) Computational results on optimal power flow show that 
the QC relaxation may bring significant benefits in 
accuracy over the SOC relaxation, especially for tight 
bounds on phase angle differences, for a reasonable loss 
in efficiency. In addition, they show that, with existing 
solvers, the SOC and QC relaxations are significantly 
faster and more reliable than the SDP relaxation. 

There are two natural frontiers for future work on these 
relaxations; One is to utilize these relaxations in power sys¬ 
tem applications that are modeled as mixed-integer nonlinear 
optimization problems, such as the Optimal Transmission 
Switching, Unit Commitment, or Transmission Network Ex¬ 
pansion Planning. Indeed, Mixed-Integer Quadratic Program¬ 
ming solvers are already being used to extend these relaxations 
to richer power system applications J4j, (44) . (45) , (46), (47) . 
The other frontier is to develop novel methods for closing 
the significant optimality gaps that remain on a variety of test 
cases considered here. 
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Appendix 

Proof of Femma IyTI 

Proof: 

Property 1 - the absolute square of power - 

l %| 2 = \Ya \ 2 {W 2 - WuWij - W U W* + | Wij\ 2 ) 
is derived using the following steps: 

Sij = Y*W U - Y*Wij 

S,,S*, = (YjWu - Y^Wj'hY.jW* - YijW*j) 

\Sij \ 2 = I Y^WuW*, - Y,, -U *ir„ - | YitfWaW* 

+ \ Y i: j\ 2 WijW*j 

\Sij \ 2 = l^l 2 ( W . 2 - WuWij - W U W* + | Witf) . 


Property 2 - the absolute square of the voltage product - 

\Wij \ 2 = w 2 - W*Z*S tJ - W„Z,,S* + IZ^lSif 2 


is derived using the following steps: 

Sij = Y*W U - Yj*Wij 
Wij = Wa - ZtjSij 

WijW* = (U'„ - ZtSiiXW* - Z, ; .S'* ! 

\Wijl 2 = w 2 - WfZf'S,, - WuZ^S* + \z t ^\s^. 


Property 3 - the absolute square of current - 

/„• = \Y,j\ 2 (W vl + W :jj - W tJ - W*) 


is derived using the following steps: 


hj = Y,(V, - V)) 

Iiilfi = Yij(Vi - V,)Y*(Vr - V*) 

\Iij\ 2 = I YijfiViV* - V t V* - V*V 3 + VjV?) 
kj = \YiiWWa - Wa - W* + W J3 ). 


Model 7 AC-OPF-W with Extensions 
variables: Sf (Vi e TV), W y (Vi, j e TV) 
minimize: ( |16a| > (34a) 

subject to: ( | 1 6b[ >— ( fl6d] >, ( | 1 6h| >— ( p~6il > 

S? - Sf - YfWa = ^ Sij Vi e N (34b) 

(■ i,j)eEuE R 


■y/-* 


, Wa 

~ Y*j ^ (i,j)eE 

(34c) 

2 J 


•y/'* 

ij 

<S>. 

1 O - 

1 Wjj 

w* 

Y*^ (i,j)eE 

T ij 

(34d) 


Extensions for Transmission System Test Cases 

In the interest of clarity, properties of AC Power Flow, and 
their relaxations, are most often presented on the purest version 
of the AC power flow equations. However, transmission system 
test cases include additional parameters such as bus shunts, 
line charging, and transformers, which complicate the AC 
power flow equations significantly. Model [7] presents the AC 
Optimal Power Flow problem (similar to Model |2j with these 
extensions. In the rest of this section, we show that the results 
of Section [V] continue to hold in this extended power flow 
model. 

The Two SOC Formulations: In this extended power flow 
formulation, the second-order cone constraint based on the 
absolute square of the voltage product 12 remains the same, 
i.e., 


\w i3 \ 2 < WhWjj (35) 

However, the constraint based on the absolute square of the 
power flow 0 is updated to include the transformer tap ratio 
as follows: 






Wu 

WA 211 


and the power loss constraint (|22d|> is updated to 


(36) 


Property 4 - voltage drop - 

Wjj = Wa - ZfjSij - ZijStj + \Za\%j 


is derived using the following steps: 


Wjj = Wa - Wa + W tj - W u + 11 * 

+ Wa - U , 7 - W* + Wjj 

Wjj = Wa - u „ + Wij - U ,, + W*j + I Zij\% 

Wjj = Wa - ZfjSij - ZijS*j + \Zij\%. 


Si. 


T Sji — Zij 




+ bfjQ(Sij) 





(37) 


We now show that, when the power flow equations P4c| >-( |34d| ) 
are present in the model, these two versions of the second- 
order cone constraints are also equivalent, i.e.. 


Sij = (Y* ~ if-) ~ Y*'^ (i,j) e E 

Sji = (y* - if) Wjj - Y* (i, j) e E 
\Wij \ 2 < WaWjj (i,j) g E 


Note that property 3 is used in the second step. 


(W-E-SOC) 
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is equivalent to 


S„ = Y* 


#) 


Wjj _ -y* Wjj 

\Tij\ 2 T* 

* 


/ b c \ W* 

Sji = (Y*- — Wjj — Y*j -jrj- 


Sij H - Sji ^ij 


-4( 


7 .. + ( lii. 

V 11 ' \ o 


\Si: 


< 


m 

) (hj) 

lij (^7 j) ^ ^' 


Tij j 2 + Wjj 

Wii 


Wii 

|Ty| 

e £ 


(*,j)eE 
(*,j) e ^ 

+ ^V^(.s'w)) 

(C-E-SOC) 


We begin by redeveloping the properties of Lemma | V. 1 1 in 
the extended model. 

The Equalities: As both models contain constraints ( |34c[ >— 
( |34d[ ), the properties arising from these equations can be 
transferred between both models. 

Proof: 

Property 1 - the absolute square of power - 


\SiA 2 



w iz w i0 

\Tij\ 2 T* 





Wa W* I Wij\ 2 \ 

\Tij\ 2 T tj \Tijl 2 ) 

(38) 


The derivation follows similarly to the one presented earlier 
and the details are left to the reader. The only delicate point 
is to observe that three separate terms in the initial expansion 
can be collected into 3(5^). 


Property 2 - the absolute square of the voltage product - 

„ W 2 - 

Wf 2 = (1 - b^(z l3 ))—^ - w*z*s l3 - w,;ZijSfj 

\ J ij\ 

+ \Z l3 \ 2 + (^j + WabtpiSij)j 

(39) 

The derivation follows similarly to the one presented earlier 
and the details are left to the reader. 


Property 3 - the absolute square of current - 

W,; 


hj ~ | Yij 



Wv 

rri* 

ij 


bfMSij) 


W*. 

T :‘ + w j3 

J 


(40) 


After observing that the extension of Ohm’s Law in this model 
is given by 


la ( Yij + i 


bf 


Vi 


YjjV, (i,j) g E, (41) 


the derivation follows similarly to the one presented earlier 
and the details are left to the reader. 


Property 4 - voltage drop - 

Wa = (1 - .)'” 2 - z^Sii - z,, Sfj 

+ \Zjj\ 2 (l l3 + + b^(S tJ )) (42) 


The proof follows similarly to the one presented earlier and 
the details are left to the reader. ■ 

With these core properties updated, we are now ready to extend 
the proof from Section [V] 

Theorem A.l. ( |C-E-SOC} is equivalent to ( |W-E-Sf)C| . 


Proof: 

The proof follows the one presented in Section [V] 


a) \W-E-SOC^_ 


C-E-SOC\ Every solution to ( [W-E-SOC} 
is a solution to ( |C-E-SOC >. Given a solution to ( |W-SOC| >, by 
equality (3), we assign lj 3 as follows: 


kj = \Yjjl 2 



w± 

rj-i* 

ij 


w* 

13 


+ w* 


33 



bjj^{Sij) (■ i,j)eE 


(43) 


This assignment satisfies the power loss constraint ( [37} by 
definition of the power. It remains to show that second-order 
cone constraint in ( |C-E-SOC} is satisfied. Using equalities (1) 
and (3), we obtain 


= I Yu 



Wu Wjj 
iTiA 2 T* 


_ uc Wjj . , 

\Tjj I 4 ij \T i3 \ 2 { 


\Sr 




Wu Wjj 

\Tio\ 2 T* 


^ii _ ,c Wjj , . 

4 ij \T-\ 2 


Wu W* \ Wt] \ 2 \ 

\Tij\ 2 Tjj \Tjj\ 2 J 

(44a) 

W u W* WuWjA 
\Tij\ 2 Tjj \T i3 \ 2 J 

(44b) 


\Sij\ 2 


< 







w* 

IJ 




Wu ( f bjj \ 2 Wu 
\Tij\ 2 2 ) 


I s tj \ 2 


< 






(44c) 

(44d) 


b) C-E-SOC ^>| W-E-SOC- Every solution to ( |C-E-SOC} 


is a solution to ( [W-E-SOC I. We show that the values of 
Wjj in ( C-E-SOC} satisfy the second-order cone constraint 
in ( |W-E-SOC} . Using equalities (2) and (4) and the fact that 
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Fig. 3. Runtime Profiles for the Two SOC Relaxations. 


Fig. 4. Runtime Profiles for the Two QC Relaxations. 


Wu = W*, since Wu is a real number, we have 


\W tj \ 2 = (1 - bfMZij)) 


w 2 

’ r or 


T - W*Z*S tJ - W„Z ij Sf j 


+ \Zij\ 2 


W : \ 2 \s,,\ 2 +(f 


6JL-V W 2 

3 ' “ + WubfMSij) 


\Tij\ 2 


(45a) 


W- 


\wn \ 2 < (i - 7 ." 2 - 

+ |Z 8J | 2 (ir„/, ; + 7^2 + u-„b:va(.s/)) (45b) 


|W«r < ( (1 - 


w,. 


r r 2 ZtjSii ~ Z 'J S ti 


\Wijl 2 < ir„u /; 

and the result follows. 


Corollary A.2. Model [5] z.v equivalent to Model 
extended AC Power Flow formulation from Model [7 


transmission system test cases from the NESTA v0.4.0 archive 
03 . Figure[3]compares the two variants of the SOC relaxation 
and Figure Hfcompares two variants on the QC relaxation. 

Both figures indicate that the two formulations are very 
similar for small test cases but, on the larger test cases (i.e., 
with more than 1000 buses), the C-QC formulation has a faster 
convergence rate, in IPOPT. This suggests that, despite its 
increased size, the C-QC formulation originally presented in 
0 is preferable from a performance standpoint and that the 
C-SOC formulation may be preferable on very large networks 
(e.g. above 9000 buses). 


+ bfMSij) (45c) 


(45d) 


6l in the 


Comparison of SOC formulations 

Section [V] proposed two equivalent formulations of the 
second-order cone constraints for power flow relaxations. 
Although both formulations define the same convex set, it is 
unclear if they have the same performance characteristics. For 
example, the current-based constraint ( |C-SOC| i has more con¬ 
straints and more variables than the voltage-product constraint 
( W-SOQ . All other aspects being equal, one would expect 
( C-SQC > to be slower than ( |W-SOC| >. This section investigates 
the performance implications of these two formulations on 
both the QC and SOC power flow relaxations. Four power 
flow relaxations are considered, W-SOC (Model [4]), C-SOC 
(Model [4] with flC-SOCD ), W-QC (Model[5]with | |W-SOQ ), and 
C-QC (Model [5}. To test the performance of these relaxations, 
each model is evaluated on 105 state-of-the-art AC-OPF 





































